Thermal pion masses in the second phase: \fj,j\ > 
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Density and thermal corrections to the mass of the pions are studied in the framework of the SU (2) 
low energy effective chiral lagrangian, in terms of the isospin chemical potential . We concentrate 
the discussion in the region where the isospin chemical potential (absolute value) becomes bigger 
than the pion mass at zero temperature and density, i.e. in the phase where the condensed n~ 
phase appears (for negative chemical potential). We are able to calculate the thermal and density 
evolution of masses in the limits where 3> m and where \fii\ > m. We also identified the phase 
transition curve. 
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In this paper we continue with the discussion of the 
thermal and density pion properties that we started in 
Q Q now in the so called second phase where the ab- 
solute value of the isospin chemical potential becomes 
bigger than the tree-level pion mass at zero temperature 
and density. In the first phase, unless we go to the chi- 
ral limit, the Isospin symmetry SU(2) is broken due to 
mass effects. This effect, as soon as we take |^/| > 
becomes extremely important since one of the charged 
pions condense. Here we will use the following definition 
of the isospin chemical potential /ij = fj, u — fid-, being [i u ^ 
the chemical potential corresponding to the baryon num- 
ber of the up and down quarks. In the case of pions, the 
isospin chemical potential will correspond essentially to 
the difference between the number of positive and nega- 
tive charged pions. 

We recall that pion masses at finite temperature m w (T) 
have been studied in a variety of frameworks, such as 
thermal QCD-Sum Rules 0, Chiral Perturbation The- 
ory (%PT) (low temperature expansion) Q], the Linear 
Sigma Model [], Q, the Mean Field Approximation J], 
the Virial Expansion ||, etc. In fact, the pion propa- 
gation at finite temperature has been calculated at two 
loops in the frame of x?T H [l(J . In our article, we follow 
the way of introducing chemical potentials in xPT that 
was presented for the first time in [jn], Q. Even though, 
these articles deal with QCD with two colors rather than 
QCD with three colors, it is clear that both problems are 
intimately related. 

The introduction of in-medium processes via isospin 
chemical potential has been studied at zero temperature 
13], 0, in both phases (|^ 7 | §5 to t ) at tree level. In 
medium properties at finite density have been discussed 
for a variety of phenomena, as for example, the chiral 
condensates 1^], the anomalous decays of pions and 
etas As the density increases both, the quark condensates 
and the decay rates diminish. 

Interesting results concerning the structure of the QCD 
phase diagram, including temperature effects have been 
achieved for the case when we have simultaneously var- 
ious baryon chemical potential and isospin chemical po- 
tential simultaneously. Two completely different ap- 
proaches have confirmed a qualitative change in the phase 
diagram as soon as the isospin chemical potential starts 



to grow [ Ls|, |l!J . These problems concerning the structure 
of the phase transition diagram, have been also handle 
in the frame of two color QCD in four dimensions | 
as well as in three dimensions The problem with 

baryonic chemical potential has been considered in the 
frame of xPT [ {22| , and also using the finite pion number 
chemical potential 1 2 



Different properties of QCD (or QCD-inspired mod- 
els) under these kind of circumstances have also been 
analyzed in the lattice approach. For QCD with three 
and two colors, an extensive work has been carried out 
24|, [25|, ^o|, ^7], in connection with the behavior of 
different order parameters in several phase transitions in 
the \xjT plane, as for example, the transition to a diquark 
phase, the chiral condensate, etc. 

A common feature of all approaches in the second 
phase is the fact that the physical pion states do not 
correspond anymore to the usual pion charged states. In 
fact as we will see, the 7r ± states will mix in a non trivial 
way. The result of this mixture produces an extremely 
cumbersome propagator, a matrix 2 x 2, so we need to de- 
fine some criteria how to handle this propagator, accord- 
ing to the value of the chemical potential, for computing 
radiative corrections. 

This paper is organized as follows: First we briefly 
present the chiral lagrangian we will use, specifying the 
lagrangians of second and fourth order, in the momentum 
expansion. Then we present our criteria to expand the 
mass corrections in a Taylor series, in two different limits: 
a) when the chemical potential is a little bit bigger than 
the the pion mass and b) in the region where the chemical 
potential is much bigger than the pion mass. In both 
limits we will work at the lowest order in our expansion. 
This will allow us to make a selection of the vertices that 
will appear in the different diagrams. We have to note, 
however, that our procedure is only valid for values of 
the chemical potential up to the mass of the r\ meson. 
Beyond this point an SU (3) chiral lagrangian is needed. 
After the selection of relevant vertices, we compute the 
thermal mass corrections from the propagators at the one 
loop level. We employ the usual momentum or energy 
expansion of the propagators around the tree level mass. 
This will allow us to comment on the phase transition 
due to the pion condensation and compare with similar 
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results from other authors, 
conclusions. 



pCfl. Finally we present our 



I. THE CHIRAL LAGRANGIAN 

The procedure we will follow is basically the same of 
[jl], ||. The only modification will be the inclusion of 
non trivial vacuum expectation values for the pion fields. 
Nevertheless, this is not just a little detail since the dis- 
cussion becomes much more technically involved and sub- 
tle. 

In the low-energy description where only pion degrees 
of freedom are relevant, the most general chiral invariant 
lagrangian at the second order, 0(p 2 ), according to an 
expansion in powers of the external momentum is given 
by 



f 2 

C 2 = —Tr [(D^D^U + 2BM{U* 



U)} (1) 



with 



(2) 



where U is the vacuum expectation value of the field 
U. M = diag{m u ,md) is the quark mass matrix and B 
in the previous equation is an arbitrary constant which 
will be fixed when the mass is identified setting (m u + 
rrid)B = m 2 , where m denotes the bare (tree-level) pion 
mass. We will use m, to denote the pion masses after 
renormalization. 

The most general 0(p ) chiral lagrangian has the form 

U = \k(Tr[(D^D^U]f 

+\l 2 Tr [(D^D V U] Tr [(D^D V U] 
+±(l 3 + l 4 )B 2 (Tr [M{W + U)]) 2 
+\UBTr [{D^D^U} Tr [M(tf + U)] 

-\l 7 B 2 (Tr [M(W -U)]) 2 

+constants. (3) 



We have used the constants k introduced by Gasser 
and Leutwyler, [|(| for an SU(2) Lagrangian. They have 
to be determined experimentally and are also tabulated 
in several articles and books. This Lagrangian includes 
the chemical potentials in the covariant derivatives and 
in the expectation value of the U matrix. The constants 
li include divergent corrections which allow to cancel the 
divergences from loops corrections 



h(A) 



with the ms pole 



A = 



32tt 2 



li - A + In ( — 



(4) 



ln47r + r'(l) + 1. 



(5) 



For \fij\ > m there is a symmetry breaking. The vac- 
uum expectation value that minimizes the potential, cal- 
culated in , is 



U 



c - 



+ i \t\ cos < 
ins, 



T2 sin ( 



1 - 



(6) 



where c = ■ s = \ 1 — . From now, we will refer 
Hi V 

with a tilde to any vector rotated in a 



angle 



V i = V\ COS <f> + V-2 Sill 

V2 = —Vi sin 4> + V2 cos ( 



(7) 



The expanded Lagrangian, keeping all the details, is 
given by 
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£2,3 = j \Hi\s[(it 2 + K 2 3 )d ir 2 " ^M/lcTfiTT 2 ] 

£2,4 = Ji\{[- ( 97r ) 2 + M?(s 2 ^ 2 + n 2 3 - |c 2 7r 2 ) + 2Mc(t^ 2 - t^^tt 2 + (tt • 9tt) 2 } 
£4,1 = y/ifc S [4s 2 (/ 1 +Z 2 )-2c 2 ; 3 -s 2 Z4]^ 1 

£4,2 - -^M 2 {2s 2 [(d7r) 2 + 2(S t ^ 2 ) 2 ]Z 1 + 2s 2 [(cU) 2 + (c^ 2 ) 2 + (d ir 2 ) 2 ]l 2 + tfc 2 el d ir 2 l 7 
-2s 2 [qii^dtidoTt, + M?(s 2 7T 2 - 3c 2 fr 2 - ^r 2 )] (h + l 2 ) + tfc 2 [s 2 n 2 - c 2 tt 2 ] l 3 
+ [c 2 (div) 2 + 2| Mj |c(1 - 3c 2 )^^ 2 - M 2 c 2 {(s 2 + 1)tt 2 - (c 2 - s 2 )^ 2 - ^ 2 2 }]Z 4 }, (8) 



where the index n,m refers to 0(p n ), 0(ir n ) and e 2 d = where 
(m„ — rrid) 2 / (m u + m^) 2 . Note that we use a negative 
chemical potential (/ij = — as it is the case in neu- 
tron stars where there exist a condensate of ir_ cuasipar- 
ticles. The same happens for the ir + changing the sign 

Of fXj. 

From this lagrangian we are able to extract the differ- 
ent vertices that will participate in the diagrams respon- 
sible for mass renormalization. Note that in the second 
phase, there will appear vertices with three legs from 
(£2,3) and one leg from (£41). The last one is respon- 
sible for the counterterms of the tadpoles, but in the 
approximation we will use, these will not be considered, 
since tadpoles, being of higher orders in the expansion 
parameter, will be absent as we will see. 



II. SELF-ENERGY IN TWO DIFFERENT 
LIMITS: \n,\ > to, AND |/^| > to 

In xPT the natural scale parameter is A x = 47r/, 
where / is a parameter which at tree level coincides with 
the pion decay constant. Usually A x is compared with 
the tree level pion mass defining in this way the small- 
ness parameter for perturbative expansions (a = (j^j) 2 ). 
Now, since > to in the second phase, we will re-define 
the smallness parameter as 



(9) 



This is possible for values of the chemical potential less 
than m v . For higher values in energy parameters we need 
to consider the SU(3) case. 

The inverse of the free propagator for pions in momen- 
tum space, i.e. the equations of motion extracted from 
£2,2 is given by the matrix 



P~ 



tfs 2 



iD 



-2i\fij\cpo 




2i\fij\cp 
P 2 




(10) 



D 



i / d 4 xe ipx 



5 2 S[it] 
Sir ; Sir i 



(x). 



(11) 



The masses are defined as the poles of the determinant 
of the propagator matrix at zero 3-momentum, i.e. the 
solution of 



\D- 1 (P)\ P = = ° 



(12) 



We can identify the ir^ = ttq field, as it is indicated 
in |l3[ |, through its mass at = to and because it is 
diagonal in the propagator matrix. Therefore, there will 
be no difficulties to handle this propagator. This is not 
the case for the charged pions which arc mixed in a non 
trivial way. 

The free propagator for the charged pions in momen- 
tum space is 



D 



-2i|^|cp 

[p 4 - p 2 tfs 2 - 4tfc 2 p 2 } V 2i|^|cpo P 2 ~ M?s 2 



(13) 

with Dij = Dij^. The Dolan-Jackiw propagators can 
be constructed with the general formula 



ir),p) - D(k - irj,p)] 
h - Po - it 
+n B (po)[D{p + ie,p) - D{p - ie,p)]. 

(14) 

The main difficulty with this matrix propagator is that 
it is very cumbersome to integrate the different loop cor- 
rections. This fact motivated us to proceed in a sys- 
tematic way, through an expansion in a new appropri- 
ate smallness parameter, namely s when > m and c 
when \fij I 3> to. We realized that a similar expansion was 

Hiposed earlier by Splittorff, Toublan and Verbaarschot 

It is not difficult to realize that the different vertices 
which appear in our lagrangian will correspond to differ- 
ent powers of c or s depending on the case. Although it 
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could appear as a trivial correction, we will keep only the 
zero order in our calculations. As we will see, this proce- 
dure is not trivial at all and provides us with interesting 
information about the behavior of the pion masses as 
function of temperature and isospin chemical potential. 

If we scale all the parameters with \fj,,\ in all struc- 
tures, we have that: 



Propagators 

D{p; T, n^rn) 



;D(p;f,c) 



(15) 



£4 constants 

k(A) 



32tt 



[k - A - ln(cA 2 



(16) 



calculated, taking into account higher order vertices and 
propagators. We will consider these in a further discus- 
sion. 

Obviously, c and s must be smaller than the intersec- 
tion point -^=. The idea is to find a region around the 
intersection point that can be excluded safely from both 
expansions. This excluded area becomes smaller when n 
(the order of the expansion) starts to grow. By demand- 
ing that c n+1 ,s n+1 < \{^) n+1 1 we achieve this condi- 
tion. Note that if n = 0, we would exclude the region 
of the chemical potential where 

Tfm 2 < M? < V8m 2 . 
We remark again that by going to higher orders in our 
expansion, the excluded region becomes smaller, so this 
is the best bound. 

Note that bare masses (tree level masses) in this phase 
arc functions of the chemical potential as was discussed 



in 



Vertices 



iMd 



fm+n-4 



(p;A,c) (17) 



Integrals 



d d k 



d d k 

(27T)« 



A 



(18) 



m 



m+ = 
rri- = 



Mi I 



3c 2 



(21) 



As the neutral pion propagator is diagonal with respect 
to the charged ones, it's propagator will be the same in 
both limits at any order: 



D„ (p;T) 



M? 



2TTn B (\p \)S(p 2 - tf). (22) 



Self- energy 

E(p ;r,^ 7 ;/,TO,A) = fi 2 vo-(p ;T, c; A) (19) 

where from now on the bar on a parameter means that 
it is scaled with \fJ-i\- It is possible then to expand the 
propagator in powers of s n for|/ij| > m and c" when 

|Mi| > m - 

After the renormalization procedure, the divergent 
terms and the scale factor A + In A that will appear in 
the loops calculation due to dimcntional regularization 
and the £4 terms will cancel. Finally the renormalized 
self-energy will get the form 

n=0 

S *Ll»m = ^"E KHP;f) + a^(p)\nc]c n . 

(20) 

As we said before we will keep only the n = terms in 
the previous expansions. This approximation is certainly 
non trivial since it allows us to explore the behavior of 
the renormalized masses precisely in the vicinity of the 
transition point and also in the region of the high chem- 
ical potential values. Further corrections could also be 



III. RENORMALIZATION AND MASSES 

The corrected propagator is 

D c = D iD'SD H 

D' 1 = D 1 + iS. (23) 

Following the usual renormalization procedure, i.e. re- 
scaling the fields 7T; = \fZ\-K\ with Zi = \ -\-0(v) we have 
that 

^^ 1 (p) 1J = ZyiD- 1 ^ +0(y?) 

= ii?- 1 (p)„-E Ji (p) 4J (24) 

with Zij = ZiZj. The value of Z is chosen in such a 
way that the corrected propagator does not have correc- 
tions proportional to p 2 . 

As an example, let consider a free propagator and the 
self energy correction of the form 

iD- 1 (p) = p 2 - 2bp - a 2 

E(p) = E° + E 1 po + SV. (25) 

This expression for E is valid in the first phase and we 
will see that it is also valid in the second phase, up to 
0(s°). 
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Choosing Z = 1 + E 2 we have that the renormalized 
self-energy will be 

E fl (p) = (E +a 2 E 2 ) + (l] 1 + 26I]>o 

= S^ + S^po. (26) 



Unfortunately, in general, the self-energy S will not 
have the shape of eq.(|25|), but it will be a complicated 
function of the external momenta. As we want to com- 
pute mass corrections, it is possible to expand the self 
energy in terms of these corrections. In the rest frame, 
where p = 0, the energy will be po = m R = m t + 5m, 
where m R is the renormalized mass, mt is the tree level 
mass and dm is the correction due to the self-energy 
terms. Then 



S(m fl ) = S(m t ) + S'(m t )<5m+ iS"(m t )(<5m) 2 + 



E°[m t ] + E 1 [mt]m R + T, 2 [mt]m R + 0(5r 



(27) 



where the mass inside the brackets is the mass around 
which the self-energy expansion was computed. 

As we said in eq. (|l2"|) , the masses are defined as the 
poles of the determinant of the corrected propagator ma- 
trix at p = 0. In the case of the neutral pion, since it 
does not mix with the other terms, it is very easy to find 
the mass correction as the solution of D R 1 (po) 00 = 0. 

As we said before, due to the complicated form of the 
self-energy, we expand the renormalized mass in the rest 
frame in powers of the tree-level mass corrections 5m. 
Since the renormalized massed will be extracted as solu- 
tions of 



D R \m R ) 



0, 



(28) 



then we only need to compute the dm corrections. This 
can be done taking an expansion of the determinant up 
to 0{5m) 2 , which leaves us with a quadratic equation for 
5m. 

After finding the 5m corrections, which are of course 
given in a power series of v, we can neglect terms of 
0(v 2 ) in the renormalized mass. Note that for finite tree 
level masses (7714 ^> t'l/Uil) it is enough to expand up to 
0{5m) ~ 0(v) 

In the case of the condensed pion, where the tree-level 
mass vanishes, some care has to be taken because now 
we do not have a natural dimensionfull quantity to refer 
to as a scale parameter. Here we will find corrections 
to the mass, only for finite temperature, in the same 
way as it occurs for the photon-mass corrections in QED 
pl|, which in our ca se will include terms of the form 




FIG. 1: Propagators O(s ) 




by 



FIG. 2: Relevant vertices at O(s ) 



IV. LIMIT \fi t \ > m 



A. Propagators and vertices (\ni\ > m) 
The propagator for the charged pions at O(s ) is given 



D 



p 2 -2i\ni\p Q 
[P 4 -^ 2 P 2 ] \2i|M/bo P 2 



-0{s 2 ) 

Pi 

(29) 

It is more convenient to work with the following combi- 
nation of fields 



(30) 



We remark that the 7r± fields do not correspond to the 
physical charged pion fields but to a combination of them. 
This fact is a consequence of the non trivial vacuum 
structure in this phase. This combination is also not triv- 
ial due to derivative terms and the inverse D'alambertian 
operator. We would like to mention that in the first phase 
(pi < m) the tt± correspond effectively to the charged 
pions and the propagator becomes anti-diagonal. Then, 
the equations of motion for charged pions are 



ib 1 -- 



\As 2 



p 



2|p 7 |cp - \li 2 s 2 



p 2 - 2\n,\cp 



D~l D-l 



D- 1 D- 



\As 2 



(31) 
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and the free propagator matrix is 

1 



D = i 





i 

p' 2 +2\^ I \p a 



p 2 — 2\fii \pa 





+ (32) 



The Dolan-Jackiw propagators for charged pions, in- 
cluding isospin chemical potential (enough for one loop 
calculations) are 



and the other propagators, D ++ and D__ are of order s 2 , 
so we will neglect them. For diagrammatic purposes, the 
D + _ propagator will be drawn as an arrow from + to — 
as appear in FIG. |l|. 



D{p;T) + _ 



D(p;T)_, 



p 2 - 2|^ 7 |p + 
+27m B (|p |) ( 5(p 2 -2|^|p ) 



0(s 2 ) The relevant vertices at the zero order in this region 

are shown in FIG. ||, where the corresponding analytical 
(33) expressions arc given by 



T/OOOG 

V 2A 



y00± 
V 2A 



v 4 



4.2 

— h 
4.2 



jp ( 2 P(o ) '?(o) +2p (+) p ( _ ) + \lM\(p° +) -P°-,) - Mi 2 - [P(o) +9(0)] +P(->] +/i?C(s 2 )} 
pa { [P(+) + [P(-) + 3c-)] - 6 M? - 4|mi|(p° +) + ~ P°-) - Q°-)) ~ 2p (+) q {+) - 2p ( _ )(?( _ 
-»2^ [ M 2 Z 3 + (2 P q + M ?)/ 4 - faldh + tfO(s 2 )] 



(34) 



r 




FIG. 3: E 00 at O(s ) 



we said, however, since we are working at zero order, this 
terms will not be considered in our analysis. 



Self-energy (|/tti| > m) 



Proceeding with the relevant vertices and propagators 
All the momenta emerge from the vertices. In principle, the loops corrections to the pion's propagator matrix are 
there will appear also other vertices: v^t, ^221 V< 2A + 1 shown in FIG. || and FIG. || 



00-- „± + + ,.±- 

} 2A ! V 2A j V 2A 



of 0(s 2 ) and v 2 ,3, «4,i of O(s). As 



Defining X = A + In A, the self-energy is 



S(p) 00 = u{-tf [|A + §/ 3 ~ 2h + 21' ~ p'} + P 2 [f A - 2f 4 + §1'] + 0(m 2 .s 2 )} 

E(p)-+ = ^{ -/i 2 [±/ 3 -2l£-8J'] -2poM[fA-2[ 4 + f/£ + f/' + 4J'] + p 2 [| A - 2f 4 + | i& + | J' ] + 0(/i 2 s 2 )} 
E(p) + _ = E(-p)_ + (35) 



with 

I' = 

J' = 

IL = 



dxy x 2 — 1 riB [ ] — 1|) + . 



Note that these definitions arc almost the same ones we 
used in our previous paper |[], ||. Here, however, the term 



d,X\/ X 2 — 1 ITlg |x — 1|) 



dx^x 2 -lx 2n 2n B {\(J.i\x). 



(36) 



FIG. 4: E_ + at O(s ) 



FIG. 5: Propagators O(c ) 



that multiplies x in the argument of the Bose-Einstein 
distribution is /ii instead of m. Following the prescription 
indicated in Section III, that the renormalized self-energy 
does not depend on p 2 we have that the renormalization 
constants are 

Z Q = 1 + j/[|A - 2f 4 + §/'] 

Z± = 1 + 4§A-2Z 4 + §!£ + §/'] (37) 

and the renormalized self-energy corrections are 
Sfloo - ^u[-\h-2I' +4:I'] 
Sft(p )-+ - v{tf[-\h+2I' Q +%J'} -8|/i z |J'p } 
£>R.(Po)+- = S fl (-j5 )-+ (38) 

plus higher corrections of order fi 2 vs 2 . We will use 
these results to extract the renormalized temperature- 
and chemical potential-dependent masses. 




T/22 T/ll T/00 

v 4,2 v i,2 V <t,2 

FIG. 6: Relevant vertices at 0(c°) 



C. Masses > m) 

To extract the m^o mass, we do not need any further 
effort since it is already diagonal in the matrix propaga- 
tor, and it's renormalized self-energy is constant in the 
external momenta, i.e., it docs not depend onpo- 



|M:|{ 



1 



21 



']}• 



(39) 



For the case of Tn n ± , the vanishing determinant of the 
charged matrix propagator, provides us with a second 
order equation in powers of p\. 



pI 



^°R- + 



4cSr 



,s 2 £<i 



R-+ 



R-+ 
0. 



R-+) 



(40) 



The solution of this equation gives the ^ masses, recog- 
nizing them respectively according to the tree-level case 
if we set the self-energy corrections to zero. We can write 
the expression for the masses in a more compact way, by 
expanding around the m + term that appears in the last 
equation. Then, the masses are given by 



m + + H Ml |[-ir 3 + ^-4J' + 0(.s 2 )] 
i| (Ul |{( S 2 + I 4-ir 3 + 2^ + 8J']) 

*v[-\h + 2/q + 8J'] + 0(^ 2 s 2 )} 1/2 .(41) 



As we indicated previously, m^- can not be expressed 
as a sum of some initial mass plus corrections, since it 
does not have a tree level mass to compare with in the 
second phase. Thermal radiative corrections, however, 
will induce a non trivial behavior for the renormalized 
m„- (T, fii). It is interesting to remark that m^- does not 
have a finite real solution for 1$ > AI' + 16J', according 
to the previous equation. This means that tt~ remains 
in the condensed state for T less than a certain critical 
value. 



V. LIMIT l/i/l >m 

A. Propagators and vertices \ ^> m) 

Proceeding in the same way as in the |^/| > m case, 
the propagator of the charged pions tt 1 and tt 2 at zero 
temperature at O(c ) is 



D 



-20(c). 

h L I 



(42) 



The 7r propagator remains the same as the > 
m case in the equation (E2n. Then, the Dolan-Jackiw 
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propagators for charged pions are 

Ai(p;T, =D 00 (p;T,\fr\) + \o(c 2 ) 

pi 

I> M (p;T, | MJ |) = + 27Tn B (\p°\)5(p 2 ) + \o{c 2 ) 

(43) 

and the propagators D 12 and l) 21 are of order c. In 
the case of a propagator with zero mass, it is necessary 



to introduce a small fictitious mass as a regulator, i.e. 
D 22 = lim rj ^ (p 2 — rj 2 ). Note that in the chiral limit 
(c = 0) the fields ir 1 and ir have the same behavior. The 
other components of the propagator matrix are of order 
c. Diagrammatically, we will denote the fl u propagator 
with a line and D 22 with a dashed line, as we can see in 
FIG. ^. The D 00 propagator remains the same. 

The relevant vertices at O(c ) are indicated in FIG. 
with the analytical expressions 



v 2 y° = ^[i+o(c 2 )] 

V2T = [2^(0)9(0, + 2p w q w + Afi 2 - (p t0 , + q w ){ Pil) + q (1) ) + fJ. 2 0(c 2 )] 

V2T = 3^ [ 2 P(o)3(o) + 2p m g m + 2tf - (p (0 , + g (0 ))(P( 3 ) + 9(2)) + M?0(c 2 )] (44) 

^2,4 2 = 3^5 [2p<D«ci) + 2p w « (a) + 2^ 2 - (p (1 , + g (1) )(p (3) + g (a) ) + M?^(c 2 )] 

^2 = i4^[(Pw3to, + + (P?o,9( 0) + + tfO(c 2 )] 

^4,2 = i4^[(Pw?ci) + + (p° + M 2 )Z 2 + K^(c 2 )] 

^2 = ^[(p (2) g (2) + 2p° )9 ° )(l 1 + h) + n 2 0(c 2 )}. 




FIG. 7: E 22 at 0(c°) 



As in the |jU : | > m case, there will appear also other 
vertices of higher order in power of c: t^lj u 2°3j 

,,122 7 ,ooi2 ,,1112 1222 12 _r /n(„\ 
u 2,3> y 2,4 ) u 2,4 > y 2,4 i ^4,2' 01 w \ c i- 




FIG. 8: E 00 (or En) at O(c ) 



B. Self-Energy » m) 

As we did for the | /x z | > m case, we use the relevant 
vertices and propagators to compute the self-energy cor- 
rections. In our approximation (C(c )), the corrections 
to the D oa and D lt are the same. This happens because 
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their vertices and propagators differ in quantities of order 
higher than c. 

The loop corrections are shown in FIG. ||, for J} u 
(which are the same as those of D 00 exchanging the single 
line with the double line), and FIG. f|for D 22 - 

Note that in this region tadpole diagrams appear, 
which are absent in the previous case. However, at the 
order c°, it turns out that these tadpoles vanish, because 



the vertex is proportional to po, and the tail of the tad- 
pole does not carry momentum. At the order c, there will 
appear non vanishing tadpole diagrams (and also new £4 
corrections). Nevertheless, for high values of the chem- 
ical potential, the leading behavior will be given by our 
approximation. 

The self energy corrections are 



£(p )oo 

£(po)n 
S(po) 2 2 



*{(p 2 o- l)[|A + 2Inc - Ih Ik + I + |/ + l« 2 T 2 ] 

+ | - 4/ + 2A 1 (p ) + 16B Cpo) + 16Bi(po) - 327riCi(p )} 

£(po)oo 

2^?{po[|A- h ~ 2[ 2 + 31nc + 2 - §/ ] + 2A 2 (p ) + 16S 2 (p ) - 32«C 2 (p )} (45) 



plus corrections of order vpqc . The functions A n , B n , C n arc dehncd as 
I 



A 2 (po) 
B (pa) 
Biipo) 
B 2 {po) 
Ci(p ) 



dx\ip^x 2 + (p — l)x] In [po( x2 ~ x) + x — it 



I dxp In [pq^ 2 — x) + 1 — ie\ 
Jo 



dxxnB(\^i\x) 



p^ - 2p x -1 + ie 



dx\/ x 2 — ln^d/z^x) 
dx\/ x 2 — lriB(\Hi\x) 



(Po - xf 



Po-1 



Ipa 



riB 



2 P o 



I'D 



pI + i 



Pq — 2p X + 1 + 26 

2l + x 

pi - 2p x + le 



C 2 ( P o) = [0(po - 2) + 6(-p - 2)]p 2 y/(p /2) 2 - ln B (\p /2\ f 



(46) 



r 



As we can see from the previous equations, it is by 
far non trivial to identify the renormalized mass, as it 
was the case in the region > m. Therefore we need 
to expand the propagator, as we explained in Section 



around the tree-level masses, identifying the term 
proportional to p\, or in this case (m t + 8m) 2 . 

For the case of Efl 00 , it is enough to evaluate E 00 with 



Po = m w o = niQ + Sitlq because D~j^ Q is diagonal with 
respect to the charged pions propagator. Then we ex- 
pand around too and rcnormalizc it, according to what 
we explained in Section III. For E^ii and £_rh we need 
to evaluate it in two values, m^i = m± + 5m±. 

The renormalization constants for the different masses 
are (evaluated at these masses) 
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Z [m ] = l + 2i/{i[A + 21nc-ii-2r 2 ] - 7 + f 7r 2 g - f - 2if 21 - 16™^n B (| Ml |)} 
Z 1 [m + ] = Z +O(vc 2 ) 

Z,[m-\ = l + |,{A + 31nc-r i -2r 2 + | + | 7 r^-f 7 r^-^5-f^-6^ 21 } 
Z 2 [m+] = 1 + l\h - 2l 2 + 3 lnc + f - ^| - 8/i 10 - 2/i 12 + 0(c 2 )} 
Z 2 [m_] = l + 2j/{|A[ 1 -2[ 2 + 31nc + 2-|^-8A' 00 } 

(47) 



where the mass inside the brackets is the mass around 
which the self-energy expansion was computed. The 
function K nm is defined as 



The self-energy corrections, associated to the different 
tree-level pion masses, as function of the corresponding 
renormalized masses are 



K = 

± 1 Tj.m. — 



dx 



(n/2). 



1 m+l ' 



(48) 





) 00 




l{[6" 


14 2 T 2 




)ll 


= S/j;(m 7r +) 00 


+ G(^ 2 i 




)ll 


= -2^ 2 


{* + 


512 fi T 6 
63 ^ 77f 






= 2z/|/i : 




2 4tt 

3 9v^3 




)22 


= 0(ft 2 


vc 2 ) 





(m n o - \fi 



(49) 



Note that the self-energy corrections actually have the 
form Ei^m^;) = S^frnj] + T^m^m^i presented in Sec- 
tion |n] 



C. Masses (\m\ > m) 



Since we have already expanded the self-energy correc- 
tions in powers of the mass corrections {8m,i) n , we can 
neglect higher terms in the determinant of the propagator 
finding then the solution for drrii from the pole condition. 
For m^o and the corrections Smo , 5m + are of order 
fi z v (neglecting terms of order vc). For m„-, because 
there is no finite tree-level mass, we need to keep in the 
expansion terms of the order v 2 ; i.e. we can consider, 
in principle, ((5m_) 2 ~ /i^ 1 / 2 . Remembering that Y*° R 
of order /i 2 v, and £^ of order ^v, we can expand the 
propagator, neglecting higher order terms. 



Summarizing the power counting is: 

Smo ~ O(hju) 
5m + ~ O(fiju) 
8m- ~ 0{^v 1/2 ) 

~ 0{^v). (50) 

The resulting expressions for the renormalized masses are 
surprisingly simple. Many terms vanish, including also 
some complex contributions. 

m n o = \fk\[l + 6vl' ] 

m v + = |^|[Vl + 3c 2 + 6^] 

TYl-ff— = 0{^vc 2 ). (51) 

Note that m w - vanishes once again in this region, i.e., 
the pion condenses again, in spite of the fact that the 
thermal corrected mass started to grow near the phase 
transition point. This behavior, however, could be just a 
fictitious result from our expansion up to order c°. 
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FIG. 10: m^-(T, n^. All the parameters are scaled with m„ 



FIG. 12: Phase diagram of the condensation point of tem- 
perature versus isospin chemical potential. The dashed line 
corresponds to the approximation made by Splittorff, Toublan 
and Veerbarschot pi. 



VI. RESULTS 



To start the discussion, we would like to mention that 
near the transition point, decreases, as in the tree- 
level approximation, and this behavior is enforced with 
temperature. In this region m n o grows with both param- 
eters (see Fig. ||). 

For , according to FIG. |o|and FIG. [ll], we see that 
a certain critical temperature must be reached before the 
7r~ is removed from the condensed state. This is precisely 
the condition that determines the phase transition curve 
(Fig. [l2]). In [Q Q,we obtained the phase transition 
diagram in the fa,T space, starting from the first phase 
by extrapolating the mass evolution m„(T, fa) up to the 
point where it vanishes. This curve has the same shape as 
the one shown in ]~~| (but for different values, confirming 
the discrepancy of about 25% the authors mentioned, 
obtained within a Nambu-Jona-Lasinio model analysis) 
and indicated previously in [[l4| p"5| . 

Starting from the first phase, according to our previ- 
ous paper Q Q , for T < , it is possible to find at first 
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FIG. 13: m w o and m n + as a function of high values of the 
isospin chemical potential at different values of temperature. 
All parameters are scaled with m 
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order in T an expression for the transition line fii (T) that 
coincides with the results given in EQ. (8.91) from |po| . 
In FIG. Ell the dashed line corresponds to this approxi- 
mation, valid at low temperature. However, as we said, 
this result gives us the transition line from the viewpoint 
of the first phase. If we start from the condensed phase, 
considering the same expansion, except that now T < \fii\ 
we find that 

^^'^'(ssxw) 173 (52) 

The reader could think from FIG.[l^ that the critical 
temperature remains constant in the second phase. How- 
ever, the temperature actually rises as function of the 
chemical potential, as we can see from the previous equa- 
tion, but very slowly, and therefore we cannot appreciate 
this behavior from the values shown in the figure. This 
growing behavior of the critical temperature, is of course 
consistent with general statements about phase transi- 
tions in the Ginzburg-Landau theory and it has been 
actually measured in the lattice for two or three color 
QCD. Nevertheless, if we consider the thermal correc- 



tions in the second phase, it happens that the condensa- 
tion phenomenon starts to disappear for a certain value 
of T > 20MeV near the transition point, remaining a 
kind of supcrfluid phase which includes the condensed as 
well as the normal phase (massive pion modes). 

In the high chemical potential region, as we said be- 
fore, the 7r~ pion condenses again. For ra T o in this re- 
gion, it increases monotonically both with temperature 
and chemical potential. m^+, as the temperature and 
the chemical potential rise, becomes asimptotically close 
to m^a as was expected (see FIG. [b|). In contrast with 
the first region, the 7r + mass grows with temperature 
and a crossover occurs somewhere in the intermediate 
region of the chemical potential for the temperature de- 
pendence, since near the phase transition point we have 
m„+(Ti, fa) > m n + (T2, fa) and this behavior changes in 
the high chemical potential region in such a way that 
m OT +(Ti,jU r ) < m 7r+ (T 2 ,^ I ). 
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